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We present a constructive derivation of a time-dependent deformation functional theory - a col- 
lective variable approach to the nonequalibrium quantum many-body problem. It is shown that 
the motion of infinitesimal fluid elements (i. e. a set of Lagrangian trajectories) in an interacting 
quantum system is governed by a closed hydrodynamics equation with the stress force being a uni- 
versal functional of the Green's deformation tensor gij. Since the Lagrangian trajectories uniquely 
determine the current density, this approach can be also viewed as a representation of the time- 
dependent current density functional theory. To derive the above theory we separate a "convective" 
and a "relative" motions of particles by reformulating the many-body problem in a comoving La- 
grangian frame. Then we prove that a properly defined many-body wave function (and thus any 
observable) in the comoving frame is a universal functional of the deformation tensor. Both the 
hydrodynamic and the Kohn-Sham formulations of the theory are presented. In the Kohn-Sham 
formulation we derive a few exact representations of the exchange-correlation potentials, and dis- 
cuss their implication for the construction of new nonadiabatic approximations. We also discuss a 
relation of the present approach to a recent continuum mechanics of the incompressible quantum 
Hall liquids. 
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I. INTRODUCTION 

Density functional theory (DFT), pioneered by Hohen- 
berg and Kohn 1 , and by Kohn and Sham 2 , is now a 
standard computational tool for studying ground state 
properties of various quantum many-body systems 3 . By 
construction, DFT, even at the exact level, yields only 
the ground state energy, and the density distribution. 
Therefore it can not be considered as a full alternative 
to the other, e. g., field theoretical many-body methods 4 . 
However, the simplicity and the computational power of 
DFT-based approaches makes them practically the only 
tool to approach realistic quantum many particle systems 
at the fully ab initio level. Moreover, from experimental 
pont of view, it is frequently quite sufficient to know just 
those quantities which are perfectly accessible by DFT 
methods. 

Conceptually DFT belongs to a class of theories of col- 
lective variables. The mapping theorems of DFT allow us 
to formally get rid of the solution of the full many-body 
Schrodinger equation, and to formulate a closed theory 
that operates with one, or a few basic collective vari- 
ables (the ground state density in the original version). 
In this respect the equilibrium DFT is very similar to 
the classical hydrostatics 5 . To find the density n(x) in 
the classical hydrostatics we need to know the equation 
of state, e. g., the dependence of pressure on the density, 
P(n), which is universal for a given substance. Similarly 
for the calculation of the ground state density in DFT it 
is enough to know a universal exchange correlation (xc) 
energy functional E xc [n}. Although the exact form of 
E xc [n] is unknown, efficient practical approximations for 
it are presently available. 

Away from the equilibrium the classical hydrostatics 



is replaced by the classical hydrodynamics. A similar 
nonequilibrium extension of DFT, a time-dependent den- 
sity functional theory (TDDFT), was founded in 1984 by 
Runge and Gross 6 . Although TDDFT is still far from 
the maturity of its static counterpart, it has already 
gained an enormous popularity in different branches of 
physics 7 . A hydrodynamic interpretation of TDDFT 
has been mentioned in the original paper by Runge and 
Gross. However, it remained practically out of use for 
more than 10 years. A comeback of hydrodynamic anal- 
ogy was caused by attempts 8 ' 9 ' 10,11 to go beyond a sim- 
ple (but formally unjustified) adiabatic local density ap- 
proximation for the time-dependent xc potential. In the 
mid-90th it has been realized that nonadiabaticity of den- 
sity functionals is strongly linked to the spatial nonlocal- 
ity. In particular, any nonadiabatic xc potential must be 
strongly nonlocal functional of the density. Otherwise the 
theory would fail to satisfy a so called harmonic potential 
theorem 12 ' 13 . An ultimate connection of nonadiabaticity 
to the spatial nonlocality is frequently referred to as an 
ultranonlocality of TDDFT. 

In the linear response regime a resolution the ultra- 
nonlocality problem was proposed by Vignale and Kohn 
(VK) 8 . They suggested to switch from TDDFT to a 
time-dependent current density functional theory (TD- 
CDFT), and to consider the current j(x, f) as a basic 
variable of the theory. VK have demonstrated that the 
xc vector potential in linearized TDCDFT can be consis- 
tently considered as a local functional of the current. In 
a subsequent paper Vignale, Ullrich, and Conti (VUC) 9 
have found an elegant hydrodynamic representation of 
the VK result - the time derivative of the VK vector 
potential can be written as a divergence of visco-elastic 
stress tensor which commonly appears in the linearized 
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classical continuum mechanics. Beyond the linear re- 
sponse VUC proposed an ad hoc extension of the nona- 
diabatic VK functional, which simply adopts the lin- 
earized visco-elastic form of the stress tensor for non- 
linear dynamics. A similar assumption has been made 
in a phenomcnological approach developed by Kurzweil 
and Baer 14,15 . 

A general resolution of the ultranonlocality problem 
in TDDFT, which has been proposed recently 16 ' 17 , also 
extensively relies on ideas and techniques borrowed from 
the classical continuum mechanics. Physically, the ul- 
tranonlocality is related to the convective motion of the 
electron fluid (in the nonadiabatic theory the particles at 
a given point of space retain the memory of their previous 
positions) 17,18 . The key idea of Refs. 16,17 was to elimi- 
nate the above source of nonlocality by reformulating the 
theory in a Lagrangian frame, i.e., in a local reference 
frame moving with the quantum fluid. Since the convec- 
tive motion in the Lagrangian frame is absent, a spatially 
local description of xc effects becomes possible. Hence 
the general resolution of the long standing ultranonlo- 
cality problem can be achieved by properly changing the 
" point of view" : while the xc potential is extremely non- 
local in the laboratory reference frame, it appears to be 
almost local from the point of view of an observer mov- 
ing with a flow. The main practical outcome of this idea 
was a rigorous derivation of a nonadiabatic local approx- 
imation for xc potential - the time-dependent local de- 
formation approximation (TDLDcfA) 17 ' 18 . A connection 
of a general nonlinear local deformation approximation 
to the phenomcnological current density functional by 
VUC 9 has been established in Ref. 19. It turns out that 
the local current density approximation by VUC corre- 
sponds to a small deformation limit of TDLDcfA. Hence 
the Lagrangian formulation of TDDFT can be considered 
as a general framework for a description of nonadiabatic 
xc memory effects. Similar to the local density approxi- 
mation (LDA) in the static DFT, the TDLDcfA can serve 
as a basic local approximation in the time-dependent the- 
ory. 

There are, however, a few important restrictions of the 
formalism developed in Refs. 16,17. In these works we 
reformulated the many-body theory and TDDFT in the 
local comoving frame by making two simplifying assump- 
tions: (i) we considered a many-body system driven by 
a scalar external potential (therefore an external mag- 
netic field was excluded from the consideration); (ii) we 
assumed that xc effects can also be described by a scalar 
xc potential. One of the aims of the present work is to 
relax both assumptions, and to formulate the theory in 
a most general form. 

Such a generalization is necessary, first of all, for the 
further refinements and extensions of the simple local de- 
formation approximation. Using the approach restricted 
to only scalar potentials 16,17 we have found that in gen- 
eral the exact stress tensor in the Lagrangian frame (i. e., 
in the space seen by a comoving observer) is a functional 
of two collective variables: a symmetric Green's defor- 



mation tensor gij, and a skew-symmetric vorticity ten- 
sor Fij. In the lowest order in spatial derivatives the 
dependence on the vorticity disappears, and we arrive 
at the local deformation approximation (TDLDcfA) 17,18 . 
If we want to go beyond the local approximation, and 
construct a gradient extension of TDLDefA 20 the depen- 
dence on must be taken into account, which seems 
to be extremely demanding technically. In the present 
paper we show that this dependence is, in a certain 
sense, trivial and can be singled out by reformulating the 
theory in terms of vector potentials (both external and 
exchange-correlation). As a result we get a theory where 
a properly defined stress tensor is, at the exact level, a 
universal functional of only one basic variable, the defor- 
mation tensor g^. As a matter of fact, this work adds 
one more member to the family of nonequilibrium DFT- 
like theories - a time-dependent deformation functional 
theory. 

Another obvious reason for reformulating the the- 
ory in terms of vector potentials is the need to de- 
scribe dynamics of many-body systems in the presence 
of external magnetic fields. Recently the ideas of the 
deformation-functional have been successfully applied to 
a phenomcnological derivation of an effective continuum 
mechanics of fractional quantum Hall liquids 21,22 . The 
present work can be considered as a formal justification 
(at the level of existence theorems) of that approach. 

The general idea of this paper is quite similar to that 
of Ref. 16,17. We formulate the many-body theory in 
the Lagrangian frame, and then use this formulation to 
derive a closed theory of a collective variable, the time- 
dependent deformation functional theory. However, the 
technique we develop here is essentially different. We 
show that it is much more transparent and economic 
to work directly with the many particle wave function, 
and to formulate the problem of quantum dynamics us- 
ing a Dirac-Frenkel variational principle. Starting with 
the quantum mechanical action in the laboratory frame 
we construct its analog in the comoving frame. Within 
this approach a theory of a collective variable emerges in 
a most simple and constructive fashion. We demonstrate 
that the time-dependent deformation functional theory 
is a natural intermediate step in solving the many-body 
problem in the Lagrangian frame. Physically this can be 
interpreted as follows. The transformation to the comov- 
ing frame corresponds to the separation of the convective 
motion of the fluid, and the motion of particles relatively 
to the convective flow. The relative^ motion is described 
by the many-body wave function \f in the Lagrangian 
space, while the convective motion is determined by a 
set of trajectories x(£, t) of infinitesimal fluid elements. 
Accordingly the complete problem splits into two natu- 
ral parts. First, one solves the many-body problem in 
the Lagrangian frame. This yields the wave function \1/ 
as a unique functional of a universal geometric charac- 
teristics of the frame - the deformation tensor g^ which 
plays a role of metric in the Lagrangian space. On the 
second step we use that solution to find the trajectories 
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x(£, t) from a closed hydrodynamics-like equation. The 
second step is, in fact, the time-dependent deformation 
functional theory, which we introduce in this paper. The 
theory is formulated both in the hydrodynamic form, and 
in a more practical Kohn-Sham form. 

The structure of the paper is the following. In Sec. II 
we illustrate the general formalism using a pedagogical 
exactly solvable example of one particle quantum dynam- 
ics. In this case the quantum problem for the relative 
motion possesses an analytic solution, and the final time- 
dependent deformation functional theory is formulated in 
an explicit form. In the concluding subsection (Sec. IIC) 
we extensively discuss the main ideas and the results of 
Sec. II. This subsection is aimed at preparing a reader 
to the most general formulation of the theory given in 
Sec. III. In Sec. IIIA the Dirac-Frenkel variational princi- 
ple is used to formulate the general many-body problem 
in the Lagrangian frame. The solution of the quantum 
problem for the relative motion is analyzed in Sec. IIIB. 
We prove the basic mapping theorem which states that 
the many-body wave function in the Lagrangian frame 
is a universal functional of the deformation tensor. This 
theorem forms a basis of the time-dependent deforma- 
tion functional theory that can also be interpreted as an 
exact quantum continuum mechanics. In Sec. IIIB the 
Keldysh-contour formalism is employed to derive a closed 
variational formulation of the theory We introduce a uni- 
versal functional W[g^] that plays a role of an effective 
"elastic" action of the exact quantum continuum mechan- 
ics. A relation of this formulation to a recent magneto- 
elasticity theory of fractional quantum Hall liquid 21,22 is 
discussed. In Sec IV we present a Kohn-Sham formula- 
tion of the theory. We introduce xc potentials both in 
the Lagrangian and in the laboratory frame. We also 
derive a few exact representations of the xc potentials 
in the laboratory frame (since just these potentials are 
of practical interest). Finally, in Sec. V we present our 
conclusions. 



II. GETTING AN IDEA: QUANTUM PARTICLE 
IN THE LAGRANGIAN FRAME 



A. Quantum mechanics in a local noninertial frame 

To illustrate main ideas, and the structure of a gen- 
eral theory developed in Sec. Ill it is instructive to con- 
sider first the simplest case of one quantum particle mov- 
ing in the presence of external vector and scalar poten- 
tials, A(x, t) and U (x, t). The system is described by the 
one particle wave function *(x, t) that satisfies the time- 
dependent Schrbdingen equation supplemented with the 
proper initial condition 



ia t *(x,t)=H*(x,t), tf(x,0) = * (x), 
where H is the usual one particle Hamiltonian 
1 



H = 



2m 



Hd x -A(x,t)] 2 + [/(x,i). 



(1) 



(2) 



For our purpose it is convenient to reformulate the 
problem of quantum dynamics using a Dirac-Frenkel vari- 
ational principle. The Schrodinger equation, Eq. (1), cor- 
responds to the condition for the extremum of the action 
S{fy*, <]>] = J f Cdt with the following Lagrangian 



£[**,*] = J 



W* d t * - 



dx, 



(3) 



where d t is a symmetrized time derivative 23 : iff* ft \P = 
(W*ftW - *ft**)/2. Using Eq. (2) we represent the La- 
grangian, Eq. (3), in the following explicit form 24 



: J dx ji 



** ft * - C/*** 



2m 



[{id^-A^Wi-id^-A^]}- (4) 



Let us transform the equation of motion, Eq. (1), or, 
equivalently, the Lagrangian of Eq. (4) to a local noniner- 
tial reference frame moving with a given velocity v(x, t). 
Formally this corresponds to a nonlinear transformation 
of coordinates x — > £, x = x(|,i), where the function 
x(£,i) is a solution to the following Cauchy problem 



dt 



v(x(£,i),t), x«,0)=€. 



(5) 



Intuitively the function x(£, t) can be viewed as a trajec- 
tory of a small element of a fluid with the velocity distri- 
bution v(x, t). Accordingly the new spatial coordinate £ 
has a meaning of the initial point of that trajectory. 

The transformation of coordinates x — > £ leads to the 
following replacements in the Lagrangian 



dx^^gd£, xi -»_0 £J) 



ft ^ ft - tfd £i (6) 



where g(£, t) = det gij is the determinant of the induced 
metric tensor 



dx k dx k 



[9i. 



ftf 

dx k dx k 



(7) 



and v(£, t) is the velocity vector transformed the moving 
frame: 



5 < (€,t) = ^u J "(x(€,t),t). 



(8) 



Substituting Eq. (6) into Eq. (4), and using the defini- 
tions of Eqs. (7) and (8) we obtain the following trans- 
formed Lagrangian 



£= I Vgdz{i**d t y 



-V Vi 



v'Ai - U 



2m 



[(*% — Ai — mvi)W][(-idti - A, - muj)*]}(9) 



Here Ai(£,t) is the external vector potential in the £- 
frame: 



~ dxi 

Mtt) = ^ r A ; -(x(€, t),t). 



(10) 
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(Raising and lowering of tensor indexes are performed ac- 
cording to the usual rules, e. g., Vi — gifv 3 .) In the deriva- 
tion of Eq. (9) we regrouped terms to obtain a physically 
expected form of Ai + mvi in the kinetic energy. 

Apparently the second term in Eq. (9) plays a role of an 
effective scalar potential in the moving reference frame. 
An important observation is that this term exactly coin- 
cides with the classical Lagrangian of a particle moving 
along the trajectory x(£, t) in the presence of the exter- 
nal fields A(x, t) and U(x,t). Indeed, using Eqs. (8) and 
(10) we find 

L c ^,t) = y^+^-[/ 

= |(x(t)) 2 + i(t)A(x(i),i) - E/(x(t),t){ll) 

where dot denotes the time derivative. Obviously the 
classical Lagrangian of Eq. (11) parametrically depends 
on the initial point, £, of the trajectory. 

The Lagrangian of Eq. (10) can be simplified by intro- 
ducing a transformed wave function \t(£,t): 

*(x(£,i),i) = <r*e iS ^<)*(£,i), (12) 

where S c i(£,t) is the classical action that is related the 
Lagrangian L cl , Eq. (11): 



S c \£,t)= / L c \tt')dt'. 



(13) 



The renormalization factor g~* in Eq. (12) accounts for 
a local change of volume induced by the nonlinear trans- 
formation of coordinates. This allows us to preserve the 
interpretation of the transformed function ty(£,t) as a 
probability density in £-space 16 ' 25 . The exponential pref- 
actor in Eq. (12) gauges out the effective scalar potential 
in Eq. (9). Inserting Eq. (12) into Eq. (9) we reduce the 
Lagrangian to the following compact form 

£[**,*] = J d£{i& d t * 

-[(id? - A)s-^*]^[H% -A^g-i^ju) 

An effective vector potential Ai(£, t) in Eq. (14) is defined 
as follows 



Mi, t) = Mt t) + vi(£, t) - %s cl (£, t). 



(15) 



The transformed Lagrangian, Eq. (14), can be also repre- 
sented in a common form of the Dirac-Frenkel functional. 
It is straightforward to check that Eq. (14) is equivalent 
(up to irrelevant total derivatives) to the functional 



de, (i6) 



where H[gij, Ai\ is the Hamiltonian in the moving frame: 
H = g~H-id e A)^(-i% - Aj)g-i. (17) 



Hence the transformed Schrodinger equation takes the 
form 

id t *(t,t) = H\ gij ,A i Mt,t), *(€,0) = * (€)- (18) 

Equations (18), (17) and (15) completely determine 
the dynamics of a quantum particle in the local nonin- 
ertial frame moving with the velocity v(x, t). The cor- 
responding Hamiltonian, Eq. (17), contains two types of 
"external" fields, a tensor field g^ , and a vector field Ai . 
The tensor field gij(i,t) plays a role of an effective met- 
ric, which produces a "geodesic" inertia force 16 . This 
force makes a free particle to move along geodesies in the 
deformed £-spacc. The effective vector potential Ai(n, t) 
is responsible for a combined action of the physical ex- 
ternal forces and the rest of inertia forces (i. e. the linear 
acceleration force, and the generalized Coriolis and cen- 
trifugal forces) 16 . 

To complete the discussion of dynamics in a general 
noninertial frame, we present two fundamental conser- 
vation laws that follow from the equation of motion, 
Eq. (18) (for a detailed derivation see Appendix). These 
are the continuity equation 

d t n + = 0, (19) 

and the local momentum balance equation 

dtlk -?(% A k - d e A t ) + nd t A k + V?V^ = 0. (20) 

In equations (19) and (20) n = \^>\ 2 is the probability 
density, and j l and are the current density and the 
stress tensor, respectively. The operator Vi in Eq. (20) 
stands for the covariant derivative in a space with metric 
gij (sec, for example Ref. 26). In particular the stress 
force, ViP^, in Eq. (20) is the covariant divergence of a 
second rank tensor: 

= -j=d?y/9% - lr yij O e g u . (21) 

In Appendix we show that in general the current density 
and the stress tensor can be expressed in terms of the 
functional derivatives of the Hamiltonian H[gij,A\ with 
respect to the vector potential and the metric tensor, 
respectively: 



j\i,t)- 



(22) 



(23) 



In the case of a single quantum particle Eqs. (22) and 
(23) lead to the following explicit representations 

i — — — — 77 

h = - — (**%* - *%**) - —A k , (24) 



p . = 

13 2m 



2m - " m 

{Kig-i^)*(Kjg-i^) + c.c. 
vj>*\|)- 



1 

hi 

9ij 



2^5 ? V S V5 
where Kj — —idy — Aj is the kinematic momentum. 



(25) 
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B. Quantum particle in a comoving frame 

The form of the equations of motion, Eqs. (18) and 
(15), is invariant under the transformation to an arbi- 
trary local reference frame defined by its velocity - the 
vector- valued function v(x, t). To specify a particular 
frame we need to supply the above system of equations 
by a local "gauge- fixing" condition 16 . One of the most 
simple and natural choices of such a gauge condition is a 
requirement of zero current density: 



j(€,*) = o. 



(26) 



This requirement specifies a comoving Lagrangian frame. 
That is a reference frame moving with the velocity 
v(x, t) = j(x, i)/n(x, t), where j and n are the current 
and the density in the laboratory frame. In this case the 
new coordinates £ become the Lagrangian coordinates, 
while the metric tensor gij, Eq. (7), acquires a mean- 
ing of the Green's deformation tensor 27 . Substituting 
Eq. (26) in the the continuity equation we observe that 
in the Lagrangian frame the density distribution is sta- 
tionary and equals to the density at t = 0, which is fixed 
by the initial conditions 



*f)^(eo)^o(o = i*o(6f 



(27) 



Similarly the local momentum balance equation, 
Eq. (20), in the Lagrangian frame simplifies as follows 



n d t Ak + VgV l P l k = 0. 



(28) 



Equation (28) reveals a physical significance of the ef- 
fective vector potential Ak- It produces a force (an ef- 
fective electric field) that exactly compensates the local 
stress force. As a result the net force exerted on every 
infinitesimal volume element vanishes, which guarantees 
vanishing current and a stationary density in every point 
of the Lagrangian £-space. 

Using Eqs. (24) and (27) we can represent the gauge 
condition of Eq. (26) in the following explicit form 



2rn 4 ? ' 



no 
m 



A k - 



(29) 



A complete set of equations of motion in the La- 
grangian frame consists of three equations. These are 
the Schodingcr equation (18), the zero-current condition 
(29), and Eq. (15) that relates the effective vector po- 
tential A to the external fields. The solution of this sys- 
tem yields the wave function <£(£, t), and the trajectory, 
x(£, t), of the Lagrangian frame for a given configuration 
of the external fields, A(x, t) and U(x,t). 

The problem of finding x(£,i) from Eq. (15) can be 
brought to a more physical form. The time derivative of 
Eq. (15) takes the form 

mdm + d t Ai - % (^v k v k + v k A k - £/) - d t A = 0. 

(30) 



The first three terms in the left hand side of Eq. (30) cor- 
respond to a combination of the external and the inertial 
forces, while the last term is precisely equal to the local 
stress force [see Eq. (28)]. Hence Eq. (30) [i. e. the time 
derivative of Eq. (15)] has a clear meaning of the force 
balance equation in the comoving frame. On the other 
hand, it can be considered as an equation of motion for 
the dynamic variable x(£, t). Indeed, using Eq. (5), and 
the explicit representations for v 1 , Eq. (8), and for Aj, 
Eq. (10), one can straightforwardly reduce Eq. (30) to 
the following form 



m- 



dx l 

~dt 



dAk _ dAA f dA k 

dx i dx k J \ dt 

dU dC dA 

dx k dx k dt 



0, (31) 



where (dAk/dt) x means the time derivative at fixed x. 
Equation (31) is exactly the Newton equation for a classi- 
cal particle moving in the presence of the external electro- 
magnetic force and the stress force [the last term in the 
left hand side of Eq. (31)]. 

Thus the complete system of equations, which deter- 
mines the quantum dynamics in the Lagrangian frame, 
can be rewritten as follows 



id t *{&t)=H\g ij ,AiWt,t), 
A k = -2^(**%* - 



(32) 
(33) 



<9f 



mx k = [x x B(x, t)] k + E k (x, t) + dtAi, (34) 

where E(x, t) and B(x, t) are the external electric and 
magnetic fields, which are defined in a usual way: 

E(x,i) = -^A(x,i)-S x C/(x,i), (35) 
B(x,i)=d x x A(x,i). (36) 

The system of Eqs. (32)-(34) should be solved with the 
initial conditions 

*(€,0) = *o(€), x(€,0)=£, *(€,0)=v (€), (37) 

where Vo(£) is the velocity distribution in the initial 
state ^o(^)- The Hamiltonian H[gij,Ai] is defined af- 
ter Eq. (17), and the metric (deformation) tensor is 
connected to the solution, x(£, t), of Eq. (34) via Eq. (7). 

Interestingly, the whole system of Eqs. (32)-(34) can 
be obtained from a single variational functional of the 
following form 

C[*,A,x] = J de{*** K $ - ^*H[ gij ,Ai}^ 

+ no(€) [y (x) 2 + xA(x, t) - tf(x, t)] }{38) 

Apparently the first two conditions for the extremum of 
Eq. (38), SC/S^* = and 5£/5A = 0, are equivalent to 
the Schrodingcr equation, Eq. (32), and the zero-current 
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constraint, Eq. (33), respectively. The third condition, 
S£/5x = 0, yields the equation 

mx k = [x x B(x,t)] fc + £7 fc (x,i) - ^<*l^l*>- ( 39 ) 

By direct calculations one can check the following iden- 
tity 

where the variational definition of the stress tensor, 
Eq.(23), and the local momentum balance equation of 
Eq. (28) have been used. Hence the last term in the 
right hand side of Eq. (40) is identical to the correct 
stress force entering the equation of motion for x(£,t), 
Eq. (34). 

The Lagrangian C, Eq. (38), describes classical dynam- 
ics of infinitesimal fluid elements, coupled to constrained 
quantum dynamics in a space with metric g^ . It is worth 
noting that the coupling is purely geometric (minimal) - 
the classical trajectories x(£, t) enter the quantum part of 
the problem only via the induced metric (£, t), Eq. (7). 

The complete set of Eqs. (32)-(34) consists of two 
parts: (i) the universal, i. e. independent of external 
fields, quantum problem defined by Eqs. (32), (33), and 
(ii) the "classical" equation of motion, Eq. (34), for the 
trajectory of the comoving frame. The universal problem 
of Eqs. (32), (33) corresponds to quantum dynamics in 
the space with a given time-dependent metric, subjected 
to a constraint of zero current density. The solution of 
this problem (provided it exists and unique) determines 
the wave function and the selfconsistent vector potential 
as universal functionals of the metric tensor, ^[s'jjK^, t), 
and A[gij](£, t). Hence the stress force, the stress tensor 
Pij , as well as any other observable are also functionals 
of the metric (deformation) tensor. Substituting the so- 
lution of the universal quantum problem into Eq. (34) 
we obtain a closed equation of motion for the trajectory 

x(€,t) 

mx k = [xxB(x,t)] t +£ t (x,f)-^ViP/[ 5ij ] (41) 

Equation (41) is easily recognized as an equation of a 
nonlinear elasticity theory in the Lagrangian formulation 
of continuum mechanics 27 . In thisjeontext the functional 
dependence of the stress tensor P- on the deformation 
tensor gij plays a role of the exact equation of state. 
This equation of state is determined from the solution 
of the universal quantum problem defined by Eqs. (32), 
(33). 

In the one particle case the universal problem of 
Eqs. (32), (33) is exactly solvable. Therefore the exact 
quantum equation of state can be found in an explicit 
form. Indeed, it is easy to see that the following wave 
function and the selfconsistent vector potential 



satisfy the system of Eqs. (32), (33) if the phase ip(£,t) 
takes the form 



9 'k \99 J 'k 9 dt' (43) 



where no(£) and <fo(£) are the density and the phase 



of the initial state: ^o(^) 



n e 



1<P0 



Substituting the 



wave function, Eq. (42), into Eq. (25) we find the stress 
tensor functional (the quantum equation of state): 



9v 
4^5 



n 



V9 



(44) 



The corresponding stress force, which enters the "elastic" 
equation of motion, Eq. (41), can be calculated by tak- 
ing either the covariant divergence of the stress tensor, 
Eq. (44), or the time derivative of the vector potential 
A, Eq. (42). The result takes the form 



n 3 k 2m ^ 



: 'k x W '^ <l ^~ (45) 



*(€,*) = M)e 



A k {£,t) = dpipfot) (42) 



It is worth mentioning that the equation of motion, 
Eq. (41), with the stress force of Eq. (45) can be inter- 
preted as the Lagrangian formulation of the one particle 
quantum fluid dynamics (see, for example, Ref. 28). 

To complete the formal consideration of the one parti- 
cle dynamics, we note that Eq. (41) with the stress tensor 
of Eq. (44) corresponds to the Euler-Lagrange equation 
for the following "elastic" Lagrangian 

£ c; [x]= J ^{n [y(x) 2 +xA(x,i)-;7(x,i)] 

+ (%<T^)^(%<T^)} (46) 

The last term in Eq. (46) plays a role of a quantum elastic 
energy. As it should be, the elastic energy depends on 
x(£,i) only via the Green's deformation tensor gij(£,t). 

C. Overview of the main results and discussion 

Let us summarize main results of the present section. 
Starting from the usual Dirac-Frenkel variational princi- 
ple, we derived a complete set of Eqs. (32)-(34), which 
describes quantum dynamics in the comoving Lagrangian 
frame. This set of equations is generated by a gener- 
alized Dirac-Frenkel functional, Eq. (38), that depends 
on three functions \I>(£,i), x(£,t), and A(£,t). Two of 
them, the wave function and the frame's trajec- 

tory x(£, t), enter the theory as dynamic variables, while 
the effective vector potential «4.(£, i) is responsible for the 
zero-current constraint. This constraint ensures that our 
reference frame is indeed Lagrangian (comoving), i. e. a 
special frame where one observes no current, and a sta- 
tionary density distribution. In fact, the force produced 
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by the vector potential A, exactly compensates the stress 
force thus providing a physical mechanism of the vanish- 
ing current density [see Eq. (28)]. 

The structure of the basic Lagrangian, Eq. (38), is ex- 
tremely simple and transparent. It describes the dynam- 
ics of two coupled subsystems. These are: (i) the classi- 
cal system of infinitesimal fluid elements labeled by their 
initial positions £, and moving along trajectories x(£, t); 
and (ii) the quantum system placed in a space with met- 
ric gtj(£, t), and subjected to the zero-current constraint. 
The two systems are coupled via Eq. (7) that identifies 
the metric gij{£, t) (entering the quantum problem) with 
the Green's deformation tensor generated by the classical 
trajectories x(£, t). 

One of the most important observations is that the 
constrained quantum problem, Eqs. (32) and (33), does 
not contain any physical external potential, but depends 
only on the metric (deformation) tensor. In that sense 
the quantum problem in the Lagrangian frame is univer- 
sal - it defines the wave function as a universal functional 
of the deformation tensor, \I/(£,t) = ^[g^]^, t). As a re- 
sult any observable in the comoving frame is also a uni- 
versal functional of g^. It is, therefore, natural to call 
this formalism a time- dependent deformation functional 
theory (TDDcfFT). The physical observable of primary 
importance is the stress tensor as it determines the 
stress force in the equation of motion for the fluid ele- 
ments. The existence of the universal functional Pij [gij] 
[i. e. the existence of the solution to the universal prob- 
lem of Eqs. (32), (33)] allows us to formulate a closed 
hydrodynamics- type equation for the only collective vari- 
able, the trajectory x(£, t). It should be noted that the 
knowledge of x(£, t) is equivalent to the knowledge to the 
time-dependent density n(x, t), and velocity v(x, t) in the 
laboratory frame. The later quantities can be recovered 
from the former one as follows 



n(x, t) = 
v(x, t) = 



Mi) 

dt 



€=€(x,*) 



€=€(x,*) 



(47) 
(48) 



where £(x, t) is the inverse of x(£,t). 

In the present section we considered the simplest, but 
still nontrivial case of one particle quantum dynamics. 
An illustrative power of this example is related to the 
possibility to exactly solve the universal quantum prob- 
lem in the Lagrangian frame, and to find the universal 
functionals, *&[gij](£, t) and Pij[gij](£,t), in the explicit 
form [see Eqs. (42)- (43) and Eq. (44), respectively]. The 
exact solution of the universal problem is no longer pos- 
sible if more than one particle is present. However, as 
we will see in the next section, both the general idea of 
TDDefFT and the formal structure of the theory remain 
basically unchanged in the most general case of an inter- 
acting quantum many- body system. 



III. TIME-DEPENDENT DEFORMATION 
FUNCTIONAL THEORY 

A. Quantum many-body theory in the comoving 
frame 

Let us consider a system of N identical particles inter- 
acting via a two-body potential V(|x— x'|). The system is 
described by the A-body wave function *(xi, . . . , xjv, t) 
that satisfies the time-dependent Schrodinger equation 
with the Hamiltonian 



H = Hn 



H, 



(49) 



where Hq corresponds to the sum of one particle contri- 
butions, and Hi n t is the interaction Hamiltonian: 



H 



N 

E 

Q = l 

2^ 

a, (3 



{-id* a - A(x a ,t)) 2 
2m 

V(\ Xa - Xp \). 



+ U(Xa,t) 



(50) 
(51) 



Following the route outlined in Sec. II, we restate the 
problem in a form of the variational principle with the 
following Lagrangian 



£[**,*] = J 



i<r d t * - 



H d Xa . (52) 



The next step is to make a transformation to the co- 
moving Lagrangian frame. Formally this corresponds 
to the transformation of coordinates x Q i— » £ q : x q = 
x(£ a , t), a = 1, . . . , N, where the function x(£, t) is the 
trajectory of a fluid element, which is defined by Eq. (5). 
In addition, we introduce a renormalized many-body 
wave function ^(£1, . . . , £jv, t) in the new frame: 

N 

*(x(€i, t), . . . , x(&v, t), t) = H g-i((; a ,t)e iS «^ 

a=l 

x *(€!,...,€*,*), (53) 

where iS c i(£, t) is the classical action defined after 
Eq. (13). Equation (53) is a direct generalization of 
Eq. (12) for the N-particle system. The rest of calcu- 
lations also straightforwardly follows the line of the pre- 
vious section. Namely, we substitute Eq. (53) into the 
Lagrangian, Eq. (52), perform the abovementioned trans- 
formation of coordinates, and successively repeat all in- 
termediate steps described in Sec. II. As a result we ar- 
rive at the following generalized Dirac-Frenkel functional 
£[\&*, A., x] that describes the dynamics of A-particle 
system in the comoving frame, 

r ~ <-. „ N 
£= [z** dt - V*H[ gij ,A}y] H d£ a 

J a=l 

+ J "o«)[y (x) 2 + xA(x,i)-[/(x,t)]d£ (54) 
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where n (£) = N J U2 =2 d€ Q |*o(€, €2, • • • , 6v)l 2 is the 
initial density distribution. 

Equation (54) has precisely the same structure as the 
one particle Lagrangian of Eq. (38). The only difference 
is that the Hamiltonian H[gij,A] in Eq. (54) corresponds 
to a system of A interacting particles in the space with 
metric gij. Accordingly H[gij,A] contains two terms, 

H[ gij , A] = H [g^ , A] + [9ij] ■ (55) 

The first term, Ho, is a sum of A one particle Hamilto- 
nians: 

N | ij 

Ho [gij , A] = ^2 -9" 1 K<*,i 2m ^ a ' j9 ° l * ' ( 56 ^ 

where K aJ = -id^ - Aj{£ a ,i), and g% = g l3 (£ a ,t). 

The second term, H mtl in Eq. (55) describes a pairwise 
interparticle interaction in the Lagrangian frame, 

#int[^]=i5>(^), (57) 

a, 13 

where l£ a £ [} is the length of geodesic that connects points 
£ Q and £p in the space with metrics g^ . 

Apparently the physical significance of the Lagrangian, 
Eq. (54), remains unchanged. It describes the classical 
system of fluid elements, which is geometrically coupled 
the constrained quantum A-body system in the space 
with metrics g^ . The Euler-Lagrange equations take the 
form 

idt*(£i,...,£ N ,t) = H\gij,AiMti,...,t N ,t) (58) 

x **(€, 6, • • • , €iv, 6, • • • , Zn, t) (59) 

mx fe = [x x B(x,t)] fc + £7 fc (x,t) - ^^V^'(60) 

where P lJ '(£, i) is the expectation value of the stress ten- 
sor operator [see Eq. (23)]. An explicit expression for 
the stress tensor in the Lagrangian frame can be found 
in 16 ' . It is worth reminding that the stress force in 
Eq. (60) can be represented in three equivalent forms, 
Eq. (40). 

The most important property of the set of Eqs. (58)- 
(60) is that the equations of constrained quantum dy- 
namics, Eqs. (58) and (59), constitute a closed univer- 
sal problem, which depends only on the initial state 
^o(£i, • ■ • ,£,n), and on the metric tensor gij(£,t). The 
solution to this problem, provided it exists and unique, 
defines the many body wave function as a universal func- 
tional of * and g^, i. e. *[*o,ff*j](€ij • ■ • ,$N,t). The 
existence of this functional is the key statement of TD- 
DcfFT. It implies the existence of the exact nonequilib- 
rium equation of state, P y = P % ^[gij], and thus allows to 



formulate a closed theory of one vector-valued collective 
variable - the trajectory function x(£, t). This theory 
can be interpreted as the exact quantum continuum me- 
chanics. 



B. The universal problem and mapping theorems 
of TDDefFT 

The complete system of Eqs. (58)-(60) is simply a refor- 
mulation of the original linear Schrodinger equation with 
the Hamiltonian (49). Therefore Eqs. (58)-(60) possess 
a unique solution, provided the external fields A and U , 
and the initial state \?o are given. It is, however, not 
obvious that given metric g^ we can solve the univer- 
sal problem of Eqs. (58), (59) independently of the third 
equation, Eq. (60). The problem of well-possedness of the 
nonlinear Cauchy problem defined by Eqs. (58), (59) for 
a given gtj($, t) is equivalent to the problem of existence 
of TDDefFT. 

In Sec. IIB we have demonstrated that in the one par- 
ticle case the universal problem of Eqs. (58), (59) [for 
A = 1 they reduce to Eqs. (32), (33)] admits an exact 
analytic solution. Hence in this particular case the uni- 
versal functional *[* ,3y] (and thus TDDefFT) docs 
indeed exist. The simplest way to argue in favor TD- 
DefFT in the general case of A > 1 is to refer to the 
Rungc-Gross theorem 6 , or to its generalizations 29 . In 
fact, just this argumentation was used in the previous 
works on "geometric" TDDFT 17 ' 18 . The present formu- 
lation of the many-body theory opens up a possibility 
to approach the problem in an alternative and, possibly, 
more natural and internally consistent way. 

As we have already mentioned, within the present for- 
malism the question of existence of TDDefFT formally 
translates to the problem of the existence and uniqueness 
of solutions to the universal quantum problem, Eqs. (58), 
(59). This problem defines a map g^ which also 

assumes a map gij 1— > A, and implies the existence of 
the generalized equation of state. Substituting the con- 
straint of Eq. (59) into Eq. (58) we observe that the later 
becomes a nonlinear Schrodinger equation (NSE) with a 
special type of cubic nonlinearity. A rigorous analysis of 
the Cauchy problem for this NSE will be presented in 
a separate publication 30 . In this paper we adopt com- 
mon in TDDFT field simplifying assumptions, 6,7,29 ' 31 
and prove the uniqueness of the solution to Eqs. (58), 
(59) for analytic and v-representable metrics gij(£, t). As 
usual we call an observable v-representable if it can be 
produced in a given physical system by applying some 
external potentials. In the present context a metric (de- 
formation) tensor g^ is v-representable if it is defined by 
Eqs. (7) and (5), where v(x, t) is a physical velocity gen- 
erated by some external potentials A(x,t) and U(x,t). 

Let us assume that the deformation tensor is analytic 
in t. Hence it possesses a Taylor expansion with a finite 
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radius of convergence, 

oo (fc) 

k=i K ' 

where <7^ = d]fgij(Q). The first term in the right hand 
side of Eq. (61) is the initial metric tensor gff — g^ (0) = 
Sij, which follows from the initial condition to Eq. (5). 
Similarly one can expand the wave function $>(t) 

OO ~ 

= + ( 62 ) 
fe=i 

Now we substitute Eq. (59) into Eq. (58), and insert the 
expansions Eqs. (61) and (62) into the resulting nonlinear 
Schrodinger equation. Collecting terms with the same 
power of t we transform Eq. (59) to the following set of 
equations for the derivatives 

*W = -iF^-V, fc = l,2,... (63) 

where denote the coefficients for the Taylor expan- 
sion of the right hand side of Eq. (58). Since Eq. (59), 
which relates A. to is local in time, the nonlinear oper- 
ator in the right hand side of Eq. (58) also locally depends 
on the wave function. This implies that the fcth Taylor 
coefficient F^ contains derivatives only withp < fc. 
Hence Eq. (63) can be schematically represented as fol- 
lows 

*«=-iF< >[$>;*o], 

*W=-iF(*-i) [ff { j ?),...,^) ; * ,. ..,*(*-!)], 

Apparently this system can be solved recursively starting 
from the first equation. The solution uniquely defines a 
map: 

{dij^^o} !-» {* (fe) }. Substituting the coefficients 
$( fe ) into the Taylor expansion, Eq. (62), we obtain the 
time-dependent wave function as a unique functional of 
the metric tensor and the initial state: ^[gij, *&o](t). Im- 
portantly, this is true only if the Taylor series of Eq. (62) 
converges, i .e., if the solution to our nonlinear problem 
exists, which can not be taken for granted in general. In 
this paper we follow the common practice and make an 
additional assumption of v-representability of the metric, 
which guarantees an a priori existence of the solution and 
thus a convergence of the Taylor series. 

The above results constitute a constructive prove of 
following uniqueness theorem: For an analytic and v- 
representable Green's deformation tensor gij, the wave 
function \I>(t) in the Lagrangian frame is a unique func- 
tional of g^ and the initial state \T/o- In other words, the 
map g^, l_ » ^{t) is unique. 



This theorem is analogous to the Rungc-Gross theorem 
in TDDFT. In particular it proves the existence of the 
exact nonequilibrium equation of state, P^ = P l ^[gij], 
and the existence of a closed theory of only one collective 
variable - the trajectory x(£,t). Since the deformation 
tensor is a unique functional of the velocity, the above 
uniqueness theorem can be also viewed as a proof of the 
velocity-to-wave function mapping that forms a basis of 
TDCDFT. 

The presented proof of g^ ^ (f ) mapping is closely 
related to constructive proofs of the mapping theorems in 
TDDFT and TDCDFT by van Lceuwen 31 and Vignale 29 , 
respectively. In fact, the formal statement of the prob- 
lem in Refs. 31 and 29 is very similar to our system of 
Eqs. (58), (59). In cither proof one solves the many-body 
Schrodinger equation, supplemented by a constraint - 
the force balance equation - that relates the potential to 
a collective variable of interest (the density in Ref. 31, 
or the current in Ref. 29). The structure of the recur- 
sive calculation of the Taylor coefficients is basically the 
same in all the proofs. However, in the present formula- 
tion of the theory the proof becomes almost trivial due 
to a complete (both in time and in space) locality of 
the constraint, Eq. (59). We note that the proof of the 
Vignale's theorem can be also essentially simplified by 
reformulating the problem in a similar local fashion 32 . It 
is also important to note that neither proof attempts to 
address a question of convergence of the resulting unique 
Taylor series for the potential and/or the wave function. 
Therefore, strictly speaking, the v-representability prob- 
lem remains unresolved in any nonlinearized version of 
TDDFT, in spite occasional statements to the contrary 
in the literature. 



C. Keldysh-contour formulation of the exact 
quantum continuum mechanics 

In the previous section we have proved that the wave 
function ^(£i, . . . ,^N,t) and the effective vector poten- 
tial «4.(£,i) are universal functionals of the Green's de- 
formation tensor gij(£,t). Hence the stress tensor Pjj 
in the Lagrangian frame is also a universal functional of 
gij. In general the stress tensor Pjj is proportional to 
the expectation value of the partial variational derivative 
5H[g i j 1 A\/8gij at fixed A. [sec Eq. (23)]. However, in 
the Lagrangian frame the current density vanishes, which 
implies the following identity j k — (SH[gij,A]/Ak) = 0. 
Thus in the Lagrangian frame Pij can be also defined via 
the total variational derivative of the Hamiltonian with 
respect to the metric 

P%«]«,*) = -^(5te]lg|||5fe]), (64) 

where H[gij] = H[gij,A.[gij]]. Equation (64) relates the 
exact nonequilibrium equation of state to the solution of 
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the universal quantum problem, Eqs. (58), (59). Sub- 
stituting this equation of state into Eq. (60) we obtain a 
formally closed equation of the exact quantum continuum 
mechanics (in the Lagrangian formalism) 



mx k = [x X B(x,t)] fe +£ fe (x 1 t)-^V i P/[ S , 



(65) 



By solving this equation with initial conditions x(£, 0) = 
£ and x(£,0) = v (£) we get a set of trajectories x(£,i) 
for a given configuration of the external fields. The 
knowledge of these trajectories allows us to uniquely de- 
termine the density, Eq. (47), and the velocity, Eq. (48), 
in the laboratory frame. 

The existence of a closed continuum mechanics defined 
by Eq. (65) is a generic fact, which follows from the 
uniqueness theorem of Sec. IIIB. In spite of an appar- 
ent "classical" form of Eq. (65) , it exactly describes the 
dynamics of a quantum many-body system. All quan- 
tum and correlation effects are encoded in the equation 
of state - the functional dependence of the stress tensor 
on the deformation tensor. 

The question we address in this section concerns a 
principal possibility to formulate the exact continuum 
mechanics in a form of a closed variational principle. 
Namely, is there exists an action functional 5*[x(^, t)) that 
generates the equation of motion, Eq. (64). In Sec. IIB 
we have explicitly constructed such a functional for the 
exactly solvable one particle case. The corresponding La- 
grangian is given by Eq. (46) that defines a simple, purely 
elastic, i. e., local in time, theory A formal reason for 
this is the following identity 



<*l^l*) = /-<*l^]|*)-^ 
dgij ')</,, tig,. 



(66) 



where E[gij] is the energy functional, which turns out to 
be local in time [sec the last term in Eq. (46)]. Therefore 
in the one particle case the stress tensor is equal to the 
variational derivative of the energy functional, exactly as 
it is in the classical elasticity theory. (We note in brackets 
that the space-nonlocality, i. e. the presence of gradients 
in E[gij], is responsible for quantum effects). 

In a many/few-particle system the identity of Eq. (66) 
does not hold. In general, for a system of N > 1 par- 
ticles the stress tensor is not a functional derivative of 
any functional. Physically this is related to a relative 
motion of particles, which produces a non-instantaneous 
response of the system to a dynamic change of the metric. 
Nonetheless a variational formulation of the theory is still 
possible if one doubles the number of degrees of freedom 
by considering the evolution along a Keldysh contour 33 . 

The existence of the map gij,^o i— > ^(t) assumes the 
existence of a unitary evolution operator U[gij](t, 0) 



|tt(t)>=tf[ff«](i,0)|*o>, 
U[ 9ij ] (t, 0) = T cxp | -i J H[ gij ] (t')dt' | 



(67) 
(68) 



where T stands for the usual chronological ordering. Us- 
ing Eqs. (67) and (68) we can rewrite the definition of 
the stress tensor, Eq. (64), as follows 



V9 °9i 3 



<5 ^C/(M)|*o>. 



(69) 



Let us introduce two different deformation tensors, g^t), 
and gfj{t), and construct the following generating func- 
tional 

W\9Tj,9& =iln(* |?7 + (0,oo)C/-(oo,0)|fo), (70) 

where U ± (t,0) = U[g^](t,0) is the evolution operator 
obtained from the solution of the universal problem, 
Eqs. (58), (59), with the metric gfj{£,t). The stress ten- 
sor, Eq. (69), is recovered by differentiating W[g~j,g^j] 
with respect to g^ , and setting g^ = g^ = gij . Formally 
the operator U + (0, oo)C/~(oo, 0) in Eq. (70) describes a 
propagation from the initial time to infinity, and then 
back to t = 0. This can be viewed as a propagation along 
a closed Keldysh contour C 33 . The contour C consists of 
two branches: the "forward" (— ) branch that goes from 
to oo, and the "back" (+) branch going from oo to the 
initial time, t = 0. Using this notion one can represent 
the generating functional W, Eq. (70), in the following 
compact form 

W\gg] = iln^olIbe-'^^faSK*)*!^), (71) 

where Tq orders times along the Keldysh contour, and 
gfj takes the values g~j, and g\- on the forward, and 
back branches, respectively. The physical stress tensor, 
Eq. (69), is given by the following functional derivative 



P ij [9ij] = 



2 5W[g\ 



V9 H 



(72) 



=s«(t) 



The notation gf^ = g^ (t) means that we set metric ten- 
sors on either branch equal to the physical deformation 
tensor gij(£,t). Let us introduce the contour trajectory 
x c(£>*) that generates the contour deformation tensor 



9ij = 7J77 7J77- 



(73) 



Using Eqs. (72) and (73) one can show that the stress 
force entering Eq. (65) is the functional derivative of the 
generating functional W with respect to x c 



~ str _ dC ~ _ 5W[g%] 



(74) 



xc=x(£,t) 



Equation (74) naturally suggests the following form of 
the Keldysh action functional 



S c [x a ] 



Jo* J 



d£n 



y(x C ) 2 +X C A(x C ,t) 

-U(xc,t)] -W\g%\. (75) 
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Indeed, the stationarity condition for the action Sc[xc], 



SScixc] 



Sx c 



(76) 



x c =x(£,t) 



recovers the correct form of the hydrodynamics equation 
of motion, Eq. (65), for the trajectory x(£, t). 

The variational formulation, Eqs. (75), (76), of the ex- 
act quantum continuum mechanics is the main result of 
the present section. From the practical point of view, 
the very existence of the action functional, Eq. (75), is 
already a very useful statement. In particular it justifies 
an application of many powerful methods of the classi- 
cal continuum mechanics to dynamics of quantum many- 
body systems. One of those methods is a construction 
of an effective elastic functionals VF^y] based on funda- 
mental symmetries of a given physical system. Recently 
this approach has been used to derive a hydrodynam- 
ics theory of strongly correlated many-body states in the 
fractional quantum Hall regime 21 . One of the basic as- 
sumptions made in Ref. 21 was the existence of the ac- 
tion functional of the form of Eq. (75). The results of the 
present section provide a rigorous justification of that ap- 
proach. In the next section we will show the variational 
formulation of TDDcfFT also offers a convenient tool for 
the definition of xc potentials in the Kohn-Sham scheme. 

A few years ago a Keldysh-contour formulation of 
TDDFT was proposed by van Leeuwen 34 ' 35 to resolve 
the causality problem of the original Rungc-Gross the- 
ory. Despite certain similarities, our construction is 
fundamentally different from that of Refs. 34,35. The 
van Leeuwen's formulation of TDDFT requires the exis- 
tence of the Keldysh-contour analog of the Runge-Gross 
mapping theorem, which has not been proved up to now. 
In contrast to that, the present approach to TDDcfFT re- 
lics only on the real-time uniqueness theorem presented 
in Sec. IIIB. 



IV. A TIME-DEPENDENT KOHN-SHAM 
CONSTRUCTION 

In general the exact stress tensor Pij[gij] as well as the 
effective energy functional VF[g^] contain both kinetic 
and interaction contributions. In some systems, such 
as strongly correlated collective quantum Hall states, or 
one-dimensional Luttinger liquids, it is natural to con- 
sider the functional Pij[gij] (or VF[<?^]) as a single entity. 
This approach was successfully employed in our recent 
studies of the fractional quantum Hall liquids and liquid 
crystals 21,22 ' 36 . However, in the most of less exotic many- 
body systems, e. g., in atoms, molecules or solids, it is 
useful to extract at least a part of the kinetic contribution 
to the universal functionals, and to consider it separately 
from the rest. The Kohn-Sham (KS) construction is a 
special tool for such a separation - it allows one to calcu- 
late exactly the noninteracting part of the kinetic stress 
functionals. 



The time-dependent KS construction in TDDcfFT can 
be introduced as follows. Let us consider a system of N 
noninteracting KS particles moving in the presence of ef- 
fective potentials, As = A + A xc and Us = U + U xc . Here 
A(x, t) and t/(x, t) are the external fields, while A xc (x, t) 
and C/ xc (x, t) are selfconsistent xc potentials that are ad- 
justed to reproduce a collective variable of interest in the 
physical interacting system. In the Lagrangian formu- 
lation of TDDcfFT the proper collective variable is the 
trajectory x(£, t). In the KS system the equation of mo- 
tion for x(£, t) takes the form 



mx = x x B s (x, t) + E s (x, t) + Ff [x] 



(77) 



where F| tr [x] is the kinetic stress force that is related to 
the kinetic stress tensor, Ps,ij[gij], of noninteracting KS 
particles: 



pstr 



(78) 



In Eq. (77) E$(x, t), and Bs(x, t) are, respectively, the 
electric, and the magnetic fields associated to the effective 
potentials As(x, t) and ?7s(x, t). The functional Ps,ij[gij] 
is obtained from the solution of the universal problem, 
Eqs. (58), (59), for a noninteracting system (i. e., with 
Hi n t[gij] = 0). To determine the xc potentials one has 
to compare Eq. (77) with the corresponding equation for 
the real interacting system, Eq. (65). Apparently they 
coincide if the force produced by the xc potentials equals 
to the difference of stress forces in the interacting and 
the noninteracting systems 



E xc , k (x,t) + [x x B xc (x,i)] fc 



no dx k 



VjPLMil 



(79) 

where I* = ^/feKU) - ^[fftf ](€»*) is thc 

xc stress tensor functional, and the xc electric and mag- 
netic fields are defined as follows 

E xc (x, t) = -dt A xc (x, t) - a x C/ xc (x, t), (80) 
B xc (x,t)=a x x A xc (x,t). (81) 

Equations (79)-(81) define the xc potentials up to a 
gauge transformation. 

In practical application it is much more convenient 
to work with the KS system (i. e., to solve the time- 
dependent KS equations) in the laboratory frame. There- 
fore we need to transform the definition of xc potentials 
from the Lagrangian frame back to the laboratory one. 
This is done simply by setting £ = £(x, t), where £(x, t) 
is the inverse of x(£, t). The result of this procedure for 
Eq. (79) takes the form 

d t A xc . k - (v x (<9 X x A xc )) fe + d k U xc = ^djP xc jk (82) 

where P XCiij (x, t) is the xc stress tensor in the laboratory 
frame, which is related to P xc ,ij [gij] (£, t) as follows 



Pxc,ij (x, t) 



fla[ff«](€(x, *),*)■ (83) 
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Equation (82) recovers the force definition of the xc 
potentials introduced in Ref. 16. The most important 
new result of the present general approach is the func- 
tional dependence on the collective variables. We have 
proved that the stress tensor in the Lagrangian frame, 
(£,t), is a unique functional of only one basic vari- 
able - the Green's deformation tensor gij. The trans- 
formed xc stress tensor P XCjij (x, t), Eq. (83), which de- 
termines xc potentials in the laboratory frame, already 
depends not only on g^, but also on the function £(x, t) 
itself. However, the dependence on £(x, t) is trivial [in the 
prefactor, and in the argument in Eq. (83)], and can be 
accounted for exactly, provided the universal functional 
Pxc,ij\gij]{£,t) is known. 



Another practically important outcome is the possibil- 
ity to define the xc force, and thus the xc potentials, via 
a functional derivative of the scalar functional 
Indeed, using Eq. (74) we can rewrite Eq. (79) in the 
following form 



E xc + x x B xc = 



1 <Wxc[ffg] 



n 5xc 



(84) 



x c =x(£,t) 



W s [g^ 



is the difference of 



where W m \g$] = W\g%\ 
W-functionals in the interacting, and noninteracting sys- 
tems. The variational definition of the xc potentials, 
Eq. (84), should be more convenient for phcnomenologi- 
cal construction (e. g. GGA-like) of various approxima- 
tions. It reduces the problem to approximating a global 
scalar functional W xc [#•?], which seems to be a simpler 
task. Importantly, the very form of Eq. (84) already 
guarantees many exact constraints. For instance, any 
functional W xc [<7^;] yields the xc force that is equal to 
the divergence of a symmetric second rank tensor, which 
automatically ensures the zero net force, and the zero net 
torque conditions. It is worth outlining that the varia- 
tional definition of the xc potentials, Eq. (84), is pos- 
sible only in the Lagrangian frame. Given a functional 
W xc [gfj\ , the xc force in the laboratory frame is obtained 
as follows. One first calculates the variational deriva- 
tive in the right hand side of Eq. (84), and then makes 
the transformation to the laboratory frame by setting 
£ = £(x, t). As a result the right hand side of Eq. (83) is 
recovered. 

The W-functional of TDDcfFT also offers a con- 
venient tool for a compact, unified, and transpar- 
ent representation of all currently known approxima- 
tions in TDDFT/TDCDFT. For example, the local VK 
approximation 8 , as well as its extension by VUC 9 , corre- 
sponds to the following quadratic functional 

W™\g%] = \j d * J c dtdt'{2^ c (t-t>)SgC(t)6gC(t') 

+ [Kg(t - f ) - |/4(t - t'^g^Sgfjit')}, (85) 

where Sgij = g^ — 8^ is the linearized strain tensor, while 
K^ c (t — t') and /i xc (t — t') are, respectively, the nona- 
diabatic xc bulk and shear moduli of the homogeneous 



system, defined on the Kcldysh contour. The VK approx- 
imation given by Eq. (85) is valid in the limit of small 

A nonlinear elastic local deformation approximation 
(LDefA) introduced in Ref. 17 is generated by a com- 
pletely local W-functional of the following form 



VF 



LDof'A 



[g?j} = j d£ JdtEnigVfot)), (86) 



where the nonadiabatic xc energy density E xc (g^ ) is de- 
fined as follows 

B« : (ff«) = E{f <J '^/xc(P;no) 
p 

+ ^ff ( \ ''"/'</'•) 22(p;n )} (87) 

In this equation V(q) is the Fourier component of the 
interaction potential, f xc (p',no) is the correlation part of 
the one particle distribution function, and g2(p;n ) is 
the Fourier component of the pair correlation function. 
Both / xc (p;^o) and c^fe^o) are calculated for a homo- 
geneous system with the density no(£). In the limit of 
small deformations, when g^ slightly deviates from 5ij, 
the nonlinear elastic approximation defined by Eqs. (86), 
(87) reduces to the high-frequency limit of VK approxi- 
mation, Eq. (85). 

For comparison we also show the functional W^f 1 ^] 
that corresponds to the adiabatic local density approxi- 
mation (ALDA): 



W. 



Ad 



(88) 



where £ xc om (n) is the usual ground state xc energy den- 
sity of the homogeneous system. 

Finally we derive one more exact representation of the 
xc force. Namely, we relate A xc and U xc to the effective 
vector potential A. that enters the universal problem of 
Eqs. (58), (59). According to the momentum balance 
equation in the Lagrangian frame, Eq. (28), [see also the 
identity of Eq. (40)], the divergence of the stress tensor 
is equal to the time derivative of the effective vector po- 
tential A~ Hence the xc force in the Lagrangian frame, 
Eq. (79), can be also represented as follows 

£ XCife (x,i) + [x x B xc (x,i)] fc - -^d t A xc 4tt), (89) 

where A. xc (£,,t) = A.(£,t) — «4s(£,i) is the difference of 
the effective vector potentials in the interacting and non- 
interacting systems with the same metric g^. The right 
hand side of Eq. (89) needs to be transformed to the lab- 
oratory frame. Let us first use the standard transforma- 
tion rule to define an effective vector potential, A^^x, t), 
in the laboratory frame 



dx j 



(90) 
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Substituting Eq. (90) in the the right hand side of 
Eq. (89) we find for the stress force 

= d t A' xck \ + v^^^ + ^-A' Kci (91) 

where we used the definition of the velocity, Eq. (5). In- 
serting the result of Eq. (91) into Eq. (89) we obtain the 
following equation for the xc potentials in the laboratory 
frame 

d t A xc - v x (<9 X x A xc ) + d x ll xc = -d t A[, c 

+ vx(d*xA! KC )-d*(vA! xc ). (92) 

This equation determines the xc potentials, A xc (x, t) and 
£/ xc (x, £), up to an arbitrary gauge transformation. One 
of possible solutions to Eq. (92) takes the form 

A xc (x,i) = -.< c (x,i), (93) 
^ xc (x,t) = -v(x,t)^ rc (x,t). (94) 

Therefore there is a particular mixed gauge in which the 
xc potentials are locally expressed in terms of the effec- 
tive vector potential «4. xc . The convenience of the exact 
representation given by Eqs. (93), (94) is that it directly 
relates the xc potentials, which enter the KS equations, to 
the solution of the universal quantum problem, Eqs. (58), 
(59). Equations (93) and (94) can be useful for the anal- 
ysis of the exact properties of KS potentials, such as sym- 
metries, scaling properties, etc. 

V. CONCLUSION 

In this paper we presented a selfcontained, constructive 
derivation of the time-dependent deformation functional 
theory (TDDcfFT). The main idea of our approach to 
the time-dependent many-body problem is a separation 
of the convective and relative motions of quantum parti- 
cles. Technically these two types of motion are separated 
by the transformation to the comoving Lagrangian ref- 
erence frame. The convective motion is described by a 
set of trajectories x(£, t) of infinitesimal fluid elements, 
where £ is the initial position (the Lagrangian coordi- 
nate) of a given element. The motion of particles rela- 
tively to the convective flow is determined by the many- 
body wave function ^> in the comoving frame. Since the 
convective motion is singled out by the above transfor- 
mation, the number of degrees of freedom entering the 
quantum many-body problem is reduced. Formally the 
dynamics of the wave function ^ is constrained by a local 
"gauge" condition of zero current density. The most im- 
portant property of this constrained quantum problem is 
that it does not contain external fields. It is completely 
determined by the fundamental geometric characteristics 
of the Lagrangian frame - the Green's deformation ten- 
sor gij that enters the equations of motion as a metric 



tensor. Hence the many-body problem for the relative 
motion appears to be universal. This problem, natu- 
rally defines the wave function as a universal functional 
of the deformation tensor, ^[ffy]. Therefore the expecta- 
tion value of any observable in the Lagrangian frame is 
also a functional of the deformation tensor. In particular 
this is true for the stress force entering the equation of 
motion for the Lagrangian trajectories x(£,i). Thus the 
trajectories and hence the whole convective motion of an 
arbitrary quantum many-body system can be found from 
a closed hydrodynamics-like theory. We call this theory 
TDDcfFT since the deformation tensor is the basic vari- 
able entering all relevant universal functional. 

The set of Lagrangian trajectories x(£, t) and the cur- 
rent density j(x, t) are in a one-to-one correspondence. 
Therefore our theory can in principle be viewed as a 
particular realization of TDCDFT. However it is also 
legitimate, and perhaps even more natural, to consider 
TDDcfFT as an independent member of the family of 
time-dependent DFT-like theories, such as TDDFT by 
Runge and Gross 6 , and TDCDFT proposed by Vignalc 
and Kohn 8 . An apparent advantage of the deformation- 
based formalism is the existence of a well founded local 
approximation for xc potentials in the KS formulation of 
the theory. In fact, TDDcfFT provides the most natu- 
ral and unified framework for interpreting all currently 
known local nonadiabatic approximations 8,9,14 ' 15,17 ' 18,20 . 
The exact representations for the xc potentials derived 
in Sec. IV should result in the further progress in con- 
structing new practical nonadiabatic functionals. 



APPENDIX A: LOCAL SYMMETRIES AND 
CONSERVATION LAWS 

In this Appendix we derive two local conservation laws 
for a TV-body system placed in the space with metric 
gij(£,t), and subjected to an external field that is gener- 
ated by the four-potential Ao(£,t), A.(£,t). The dynam- 
ics of the V-body wave function . . . , £/v, t) is gov- 
erned by the time-dependent Schrodingcr equation with 
the following Hamiltonian 

N ij 

a=l 

N 1 

a—l a,P 

where g% — g % ' J '(£ a ,t) and is the length of geodesic 
that connects points £ Q and £p. 

Below we derive local balance equations from local 
symmetries of the Dirac-Frenkcl action functional (for 
a similar derivation see also Ref. 37): 

r N <-> 

S[*,Ao,A, gij \ = dt]] d{ a **(i d t -ff)*. (A2) 

J a=l 
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In particular the local gauge invariance of S is respon- 
sible for the local conservation of the number of parti- 
cles, while the general coordinate invariance of the action 
yields the local momentum balance equation. 

Let us consider the gauge invariance first. By a direct 
substitution we find that the following transformation 

$' = fe < S.«..t) ) (A3) 
A' = Ao-d t <t>(Z,t), A' l = A l + d e 4>(^t) (A4) 
S««,t)=S« (€,*), (A5) 

where </>(£, t) is an arbitrary function, leaves the action 
S unchanged, i. e., 

A' , A',dj\ = A , A, gij }. (A6) 

Inserting the infinitesimal version (<j) — » 0) of Eqs. (A3)- 
(A5) into Eq. (A6) we get the following condition 



where £/ a = £' Let us consider an infinitesimal 
transformation of coordinates, which is generated by the 
function £' = £ + rj(£,t), rj — > 0. The corresponding 
change of the action (at the extremal) takes the form 



dtd£, 



( SS . SS . . SS 



KSgij"™ ' SAo"""' ' 5 A 



SA + TT 5A l - (A15) 



In equation (A15) Sgij, 6Ao, and SAi are given by the 
infinitesimal version of Eqs. (A12)-(A14): 

Sgij = -rfdpgij ~ gikd i3 ri k - g jk dprj k , (A16) 
SAi = -r) k d e Ai - A k d e rj k - mg ik d t r) k 1 (A17) 
5A' = -r] k d e Ao + Akd t r] k . (A18) 

Inserting Eqs. (A16)-(A18) into Eq. (A15), and integrat- 
ing by parts we obtain at the following condition of the 
general coordinate invariance of the action: 



SS\ e ^ = [dtdtf-^dtt+^-deA =0 (A7) / SS SS \ . o ( r) SS SS \ 

J V SAo SA J ft ^__^_j + ^ 2ffw _+^_j 



As usual we take the variation of the action at the 
extremal "trajectory" that is defined by the equation 
SS/Sfy = 0. Therefore the change of the wave function, 
Eq. (A3), does not contribute to SS, Eq. (A7). The re- 
quirement that Eq. (A7) is fulfilled for any (f> yields the 
continuity equation 



SS 



dekAo 



SS SS . 

o^kAi — - — o^gij = 0. 



5A a ? '*"' u SA, ' Sg. 



(A19) 



Equation (A19) can be further simplified as follows 
mdtjk- j l (d^Ak - %A) + nd t A k + nd^kAo 



d t n + %f = 0, 



(A8) 



+ % 2g kj 



SS 
' Sg tj 



SS 

Sgij 



d (k g tJ = 0, (A20) 



where the density n(£,t) and the current j 4 are de- 
fined as follows 



SS SH 

n = -SA=WsAo l * h 
SS SH 

J= SAi = -^SAi^ 



(A9) 
(A10) 



where we have used the continuity equation, Eq. (A8), 
and the definitions of the density, and the current, 
Eqs. (A9) and (A10). The last two terms in the left hand 
side of Eq. (A20) are easily recognized as a covariant di- 
vergence the following symmetric second rank tensor 



Similarly the momentum balance equation follows from 
the invariance of the action S, Eq. (A2), under a general 
nonsingular transformation of coordinates, = 
where is an arbitrary (invertible) function. The 

transformation of fields, which leaves the action invari- 
ant, is the following 



2 SS 



2 ,Ti^m 



•4>(€',t) = 



N 

n 

a=l 

d£ k d£ p 



V9 Sgij Sgij 
Indeed, using the definition of Eq. (A21) we find 
SS \ SS 

1 



(A21) 



% {^Sg-J-Sg-^ 9 - 
( 1 



nm%=i,t), (ah) 

g kp (Z,t), (A12) 
^(€,*)-m^,«',t)^(A13) 



dt 



-Ao(U)+A^\t) 
2 3i3^ > z > dt dt 



dt 



(A14) 



Therefore the condition of the general coordinate invari- 
ance, Eq. (20), takes the standard form of the local mo- 
mentum balance equation 

mdtjk- f (d? Ak - % A) + nd t A k 

+ nd 6 kA a + JgViPl = 0. (A22) 

Apparently the tensor P y defined after Eq. (A21) plays 
a role of the physical stress tensor. 
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